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Upper critical field H c2 calculations for the high critical temperature superconductors 

considering inhomogeneities 
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We perform calculations to obtain the H C 2 curve of high temperature superconductors (HTSC). 
We consider explicitly the fact that the HTSC possess intrinsic inhomogeneities by taking into 
account a non uniform charge density p(r). The transition to a coherent superconducting phase at 
a critical temperature T c corresponds to a percolation threshold among different superconducting 
regions, each one characterized by a given T c (p(r)). Within this model we calculate the upper 
critical field H C 2 by means of an average linearized Ginzburg-Landau (GL) equation modified to 
take into account the distribution of local superconducting temperatures T c (p(r)). This approach 
explains some of the anomalies associated with H C 2 and why several properties like the Meissner 
and Nernst effects are detected at temperatures much higher than T c . 

PACS numbers: 74.72.-h, 74.20.-z, 74.81.-g 



I. INTRODUCTION 

High critical temperature superconductors (HTSC) 
have been discovered fifteen years ago pj, but many of 
their properties remains not well explained. Some of their 
features are completely different from the low tempera- 
ture superconductors, for example, the H-T phase di- 
agram of the HTSC possess, in certain cases, positive 
curvature for H C 2(T), with no evidence of saturation at 
low temperatures El tl These lack of saturation at 
low temperatures may minimize the importance of strong 
fluctuations of the order parameter which was one of the 
earlier leading ideas 0, Furthermore, this behavior 
cannot be explained by the classical WHH (Werthamer- 
Helfand-Hohenberg) theory for superconductors 0, Q 
where H C 2(T) exhibits negative curvature and saturates 
at low temperatures. 

In order to interpret this unusual behavior and fur- 
nishes ground to the H C 2 positive curvature many dif- 
ferent possibilities were proposed like, for instance: the 
interplay between Anderson localization and supercon- 
ductivity which affects the coherent length near the mo- 
bility edge 0, the effect of correlations among magnetic 
impurities which leads to a weakening of the pair break- 
ing process mesoscopic fluctuations [ll|, depletion 
of vortex viscosity or local pairs confined to small su- 
perconducting regions coupled through a normal host 
due to anisotropic quasiparticle relaxation rate around 
the Fermi surface |l2l. effective two component (boson- 
fermion) model 13], Luttinger liquid behavior in the 
normal state |l4l Il5j|. bipolaron superconductivity [lr^ . 
charge-density waves and stripes |l7j . 

On the other hand, another feature that has recently 
received considerable attention is the intrinsic inhomoge- 
neous distribution of charge (or doping level) and super- 
conductin g gaps, as demonstrated by the STM/S experi- 
ments |T5llf3Ll2f| . Based on these measurements and in a 
number of others, like neutron diffraction 



we have proposed a general theory to explain the main 
features of the HTSC phase diagram 25, 26]. The main 
point is that near the pseudogap temperature T* |27j 
some localized regions become superconducting, and the 
size of these regions increases as the temperature is low- 
ered and, at the critical temperature T c , these supercon- 
ducting regions percolate, and the system becomes able to 
hold a dissipationless current. This scenario is consistent 
with the Meissner effect above T c measured by Iguchi et 
al. 123 and the diamagnetic signal above T c which has 
also been measured by other groups by susceptibility [2j| • 
The drift of magnetic vortices expected to occur in the a 
superconductor phase has also been measured in HTSC 
at temperatures above T c through the Nernst effect [3(| • 

In this paper we develop a unified view for all these 
anomalous properties to explain the observed features 
of the upper critical field H c2 unusual properties in 
HTSC. In our approach the pseudogap temperature T* 
is where small static superconducting regions starts to 
develop [2H . The intrinsic inhomogeneities in the charge 
produces superconducting gaps which varies locally in- 
side a given sample and, as the temperature decreases 
below T* , the superconducting regions start to grow as 
it has been demonstrated by the Meissner effect mea- 
surements 28]. At temperatures below T*, but above 
the critical temperature T c , there is no long range order 
but there are several isolated local superconducting re- 
gions inside the sample. This scenario explains why there 
is some indications that H C 2 does not vanish at T c [sil ] 
but at a much larger value; also a much larger magne- 
tization than the expected, above T c , has been recently 
measured |29j . The main difficulty to perform calcula- 
tions in this scenario is the lack of information on the 
inhomogeneous distribution of charge in a given sample. 
Based on the experimental results which led to the idea 
of stripes [321 133| and on the STM/S results, we have 
proposed a bimodal distribution for the charge inside a 
given HTSC |26(. Although we do not know the exact 
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form of such distribution, it mimics the antiferromag- 
netic (AF) insulator regions (insulator branch) and the 
metallic regions (metallic branch) and it contains several 
features found in the cupratcs. Here wc will apply such 
distribution derived in Ref. [26j to calculate the contribu- 
tion of resulting superconducting regions to estimate the 
magnetic response for some materials and compare with 
the experimental results. This is accomplished through a 
simple model to compute the contribution of the super- 
conducting regions to the upper critical field H c2 , consid- 
ering that each of these regions acts as independent su- 
perconducting regions whose H C 2 is described by the GL 
equation. All the contributions from different regions are 
computed to give the total upper critical field of a given 
superconducting compound. This is a clear phenomeno- 
logical approach that, according to the GL theory, yield 
good results near the critical temperatures. 

It is worthwhile to mention that Ref. |l2j has also in- 
ferred that the positive curvature of H C 2 could be due to 
pair formation in small grains with local T c higher than 
the bulk T c . Therefore, we can see that several attempts 
had been published [T3, EH E2j EH connecting H c2 and 
its features with the effects of some type of disorder. By 
the same token, we apply below our recently developed 
theory [2(j on the inhomogeneities of cuprates and pseu- 
dogap phenomenon to evaluate a theory for the H C 2 of 
these materials. 

This paper is divided as follows: In Sec. II we present 
the density of charge distribution and the phase diagram 
of a selected compound of the LSCO family. In Sec. Ill the 
upper critical field H C 2 from the GL theory is considered 
and generalized with the inclusion of the inhomogeneous 
superconducting regions. In Sec. IV we compare the the- 
oretical results with some selected experimental data of 
the LSCO family and optimum Bi2212 high-T c . A good 
qualitative agreement is observed. In Sec.V we make the 
final analyses and conclusions. 



II. THE DENSITY OF CHARGE 
DISTRIBUTION 

Just for completeness, we will briefly outline the basic 
ideas concerned with the inhomogeneous charge distribu- 
tion introduced in Ref. 

To model this inhomogeneous medium we consider a 
phenomenological distribution of probability P(p(r)) of 
a given local charge density p(r). The differences in the 
local charge densities yield insulator and metallic regions. 
For simplicity we hereafter make p(r)=p. 

The distribution P(p(r)) we consider is a combination 
of a Gaussian and a Poisson distribution, which becomes 
the appropriate distribution to deal with the high and 
low density compounds 26], that is, the whole phase di- 
agram. The main features of P(p) is that it has two 
branches: an insulating one with < p < 0.05, and a 
metallic one which starts at p m . 

For most compounds p m ~{p), where (p) is the aver- 
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FIG. 1: a) The distribution P(p) for the LSCO (p)=0.15 com- 
pound is shown. The arrows indicates the density p c =0.05, 
the percolation density p p =0.225, and p m , with p m ={p). 
The hachured part indicates the region below the percolation 
threshold p p , i.e. the region below the superconducting crit- 
ical temperature T c . b) The theoretical local superconduct- 
ing critical temperatures T c (p) for this compound is shown 
as the solid line, together with the experimental pseudogap 
temperatures T* of Refs. [34l l35|l . The filled black circles are 
the experimental critical temperatures T c of the LSCO fam- 
ily from Ref. |3q|. The vertical and horizontal lines indicates 
T c (p p )=T c . 



age density of a compound. In Fig. ^ we show the 
phase diagram for the LSCO (p)=0.15 compound to- 
gether with the charge distribution and experimental 
data of Ref. 0, |3^| . Since the local critical tempera- 
ture T c (p) is a decreasing function of p, the maximum T c 
is T c ((p)), which is therefore the system pseudogap tem- 
perature T* . Upon cooling below T* part of the metallic 
regions become superconducting. At T=T C 60% of the 
material is in the superconducting phase and we say that 
the superconducting regions percolate [26j. 

In order to study the effect of the charge distribution 
on our results we have considered a constant and a lin- 
ear decreasing distribution, with metallic and insulating 
branches, evaluated with the densities p c and p p (see 
Fig^l for p m =0.15, and in the same range. The results 
are shown in Fig|2i of section IV and indicate that the 
qualitative behavior remains the same, i.e., the inhomo- 
geneities seem to be the cause of the positive curvature, 
but the quantitative agreement is worst than the results 
from the distribution derived by the STM experiments 
and used in our calculations. 
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III. THE CALCULATIONS 

It is well known that most of the HTSC are type-II 
superconductors [3g. For these type of superconductors 
there are two critical fields in the H-T phase diagram: 
the lower H c i and the upper H C 2. Above H C 2 the ma- 
terial returns to the normal metal state. By definition, 
one expects the superconductivity to disappear above the 
upper critical field H C 2- 

In the case of an external magnetic field parallel to the 
c-direction, i.e. perpendicular to the CuC>2 planes (ab- 
direction), the GL upper critical field is given by [37ll38t 

M 



H c2 {T) 



(i) 



where $o = hc/2e is the flux quantum and £, a b(T) is the 
GL temperature dependent coherence length in the ab 
plane [a USES- Therefore, H c2 is determined by the 
coherence length ^ a b{T) of the supercondutor, which is 
treated as a phenomenological parameter. In terms of 
the GL parameters the coherence length is given 



e ab (T) 



2m ab a(T) 



= &(0) 



T r - T 



(T < T c ) (2) 



where ^ lb (Q)=h 2 /2m a boT c is the extrapolated coherence 
length, m a (, is the part of the mass tensor for the ab 
plane and a is a constant j3(|. Using the BCS formula 
£ ~ VF/ksTc we expect a shorter coherence length for 
HTSC relative to the low temperature supercondutors 
due to their 10 times higher T c 's. However, due to the 
lower density of carriers |39| . vf in these materials is 
also small, which results in very short coherence length, 
£ ~ 10 A. A typical value for the extrapolated coherence 
length which we use in our calculations is £ a b(0) ~ 15 A 
and £ c (0) ~ 4l for the YBCO Jl and £ a6 (0) - 32A 
and f c (0) ~ 7 A for LSCO jM ||J where f c (0) is the 
coherence length in the c direction. One should note 
that £ c (0) is smaller than £ a fc(0) and is of the order of the 
spacing between adjacent conducting CuC>2 planes. 
Therefore, the GL upper critical field may be written 

as 



H c2 {T) = 



2< b (0) 



T c -T 



(T < T c ) (3) 



Let us now apply this expression to a HTSC with in- 
trinsic inhomogeneities in the charge distribution. 

When considering the inhomogeneity of the HTSC at 
temperatures below T* , isolated superconducting regions 
may exist in the form of se para ted islands even in mag- 
netic fields H > H C 2(T) |42ll43| . For temperatures above 
T c there is no long range order but, due to the various 
different local values of T c (p(r)) superconducting regions 
may exist in the form of separated regions. Here we cal- 
culate the upper critical field H C 2 for a given sample as- 
suming that each isolated or connected superconducting 



region displays a local H l c2 which is given by the lin- 
earized GL equation with effective mass tensor [33, E3 ■ 
Since a given local superconducting region "i" has a local 
temperature T c (i) with a probability Pi and local coher- 
ence length £j, it will contribute to the upper critical field 
with a local linear upper critical field H l c2 (T) near T c {i). 
Therefore, the total contribution of the local supercon- 
ducting regions to the upper critical field is the sum of 
all the i?* 2 (T)'s. Thus, applying Eq.@, the H c2 for an 
entire sample is 



H c2 {T) 



$ 



1 



iV 



2< 6 (0) W ^ 



T c (i) - T 
Tc{i) 



1 N 

= w ? (T < T c (i) < T c ) (4) 

i— 1 

where N is the number of superconducting regions, or su- 
perconducting islands each with its local T c (i)< T c and 
W = J2iLi Pi i s the sum of all the Pj's. As we have 
already mentioned, at temperatures above T c there are 
isolated superconducting regions, while below T c these 
regions percolate and the system may hold a dissipation- 
less current. Since H C 2 is experimentally measured at 
T < T C (H=Q), it is the field which destroys the super- 
conducting clusters with T < T c (i) < T c , leading the 
system without percolation. The first regions which are 
broken are the weakest ones, which have critical tempera- 
tures T c (i)'s lower than T c (H=0). The mechanism is the 
following: at a temperature T < T c most of the system is 
superconducting and a small applied field destroys first 
the superconducting regions at lower T c (i)'s, without loss 
of long range order. Increasing the applied field causes 
more regions to become normal and eventually when the 
regions with T c (i) w T c turn to the normal phase, the 
system is about to have a nonvanishing resistivity. This 
value of the applied field is taken as the H C 2 in our theory, 
and it is the physical meaning of Eq.Q. Thus, at a given 
temperature T, we sum the superconducting regions with 
T < T c (i) < T c , with their respective probabilities. 

It is important to notice that the STM experiments 
have demonstrated that the size of a region with constant 
superconducting gap is of the order of 20A 0, 0, |2(| ■ 
These results have also been obtained by /iSR experi- 
ments |45|. According to the discussion in the intro- 
duction, at T* some small superconducting isolated is- 
lands (or droplets) start to appear through the sys- 
tem. Therefore, near T* the GL approach should not be 
valid because the size of an island coherence lenght with 
T c (i)«T* is of the same order of the superconducting is- 
lands. As the temperature decreases the superconducting 
islands start to unite forming larger inhomogeneous su- 
perconducting regions. At T c they percolate, occupying 
60% of the system. At T < T c the occupied supercon- 
ducting volume is comparable to the size of the system 
and clearly much larger than the typical coherence length 
£j of the different inhomogeneous regions. Consequently, 
we may use the GL theory to calculate the upper criti- 
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cal field of these different superconducting regions which 
form the whole condensate, and the total sample H C 2 is 
the sum of these individual inhomogeneous contributions 
as in Eq.lQJ. 



Ref. [3l| of the LSCO family and near optimum Bi2212. 
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FIG. 2: Theoretical results of H c2 (T) (solid lines) of the LSCO 
series considering the distribution of Ref. |26j| together with 
the experimental data of Ref. @- The dashed line is a GL 
fitting of Eq. . In a) the results for a constant (dot-dashed 
line) and a linear (dotted line) distribution are also shown. 



IV. COMPARISON WITH EXPERIMENTAL 
DATA 

In this section we compare the results of Eq.QJ with 
the experimental upper critical field data of Ref. |2| and 
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FIG. 3: Theoretical results of H C 2(T) (dot-dashed line) for the 
(p) — 0.20 of the LSCO series together with the Nersnt signal 
measurements curve of Ref. [3~H (solid line) . The dashed line 
is a GL fitting. 



The experimental upper critical field H c2 of the HTSC 
may be obtained from the resistivity measurements as it 
is the field relative to a fraction of the "normal-state" 
resistivity [2J, |39j. The correct fraction which leads to 
the upper critical field is still controversial and there is 
no conclusion in whether H C 2 may correspond to the be- 
ginning, the middle, or the top (end) of the resistivity 
curves 0, Q . Another way to obtain H c2 is from the field 
dependence of the transport line-entropy derived from 
the Nernst signal j^l). Despite of this difficulty we at- 
tempt here to compare the theoretical results with the 
available experimental data. For this purpose we iden- 
tified some applied magnetic fields of Ref. Q and H* 
of Ref. as the upper critical fields for the selected 
compounds studied due to the reasons below. 

By definition, H onS et from the resistive measurements 
of Ref. Q is defined as the magnetic field at which the 
resistivity p first is detected to deviate from the zero in 
the p vs H curves, and this is the assumption used in 
Eq.(0} and, therefore, it is our definition of H C 2(T). Fur- 
thermore, it is reasonable to take T c from the H-T phase 
diagram as the temperatures where H onset =0. However, 
it is experimentally observed 0] that for some compounds 
Honset vanishes at temperatures much smaller than the 
known values of T c . This is because the superconducting 
transition is sharp at low fields, but spreads itself over a 
large interval of temperatures as the field increases. As a 
consequence one defines the field Hx, x being a percent- 
age of the normal state resistivity. In some cases the dif- 
ference between H onset and H90 (90% of the normal-sate 
resistivity) may be about 50% of T c 0. When this is the 
case the correct H C 2 may be between these two fields. In 
our calculations from Eq(0J H C 2 is calculated for T <T C 
and always vanishes at T c , with H c2 being in principle 
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FIG. 4: Theoretical results of H C 2(T) (solid line) of the near 
optimum Bi2212 together with the experimental points of 
Ref. 0. A GL fitting using Eq.JHJ is also shown (dashed 
line) for comparison. 



Ref. [IJ. For the near-optimum Bi2212 we compared our 
results with if 50 (Fig. of Ref. [2j, which vanishes at 
T c ps 80A' and is in accordance with the phase diagram 
of Ref. j2|J , while H onset vanishes at T « 50K. Also, 
in Fig. 0: and Fig. [3] the experimental points of H onset 
are shown for comparison. For the LSCO series a co- 
herence length of £ Q i,(0)=30A was adopted, which is in 
accordance with the measurements of Refs. |39l l4ll l46j. 
This value of £ afc (0) leads to ff c2 (0)=$/27r^ 6 (0)=32T. 
For the Bi2212 a coherence length of £ a b(0)=27A was 
considered, which is in accordance with Ref. |4fij . Sim- 
ilarly one gets ii C 2(0)=45T. These discrepancies for the 
low temperature values of H c2 (T) is clearly due to the GL 
expressions (Eqs.JTJ-GJ), wn i° n should not be used far 
from T c . This is the reason why we stop our calculations 
at temperatures below T c /3. At very low temperatures 
we do not know how to estimate the contributions of the 
islands with T c (i)^T c . 



V. CONCLUSION 



the H onset field. By the same token, for the Nernst sig- 
nal measurements of Ref. [2jJ , H* may be considered the 
upper critical field since it represents an intrinsic field 
which controls the onset of the flux-flow dissipation and 
vanishes at a temperature close to T c . Therefore, H* 
may be compared with H onset . In agreement with Wang 
et al. ,3l|, attempts to find H C 2 using the resistivity in- 
variably turn up a curve akin to H*. 

Now, in order to compare with the experimental fields 
of Ref. 0, we plot H C 2(T) with the measured H onset for 
the cases of (p)=0.15 (Fig. Jgfc) and (p)=0.17 (Fig. |3>) of 
the LSCO series. Also, in Fig. we plot the results of 
a constant and a linear charge distribution together with 
the bimodal distribution of Ref. [2^ . As one can see, 
the distributions yield very similar results, which shows 
that the calculations do not depend on the details of the 
charge distribution, although without a distribution of 
T c (?)'s we simply obtain a GL linear behavior. 

For the case of (p)=0.08 (Fig. we compared our 
results with H90 since this field vanishes at T c as 24K, 
which is the value of T c obtained from the phase dia- 
gram of Ref |26(, while H onset vanishes at T w 12iG 
In Fig. |3 one can see the results for ( J o)=0.20 com- 
pared with H* from the Nernst-signal measurements of 



We have calculated the upper critical field H C 2(T) for 
a disordered superconductor characterized by a distribu- 
tion of different local critical temperatures T c (i) at dif- 
ferent domains. We have applied a simple GL expression 
to each of these superconducting regions. With this pro- 
cedure we can explain the magnetic signals below and 
above T c . We have been able to fit the H C 2 curves de- 
rived by two different experimental procedures, namely 
the resistive magnetic fields and the Nernst signal |30| . 
In all the cases the curves exhibit a positive curvature 
which is different from the magnetically determined H C 2 
lines |4l| from pure GL (Eq.@). This positive curva- 
ture reflects the GL behavior of each individual domain 
in a disordered superconductor. Furthermore, taking the 
inhomogenous charge distribution into account several 
properties like the Meissner and Nernst effects, which 
are seen at temperatures much higher than T c , are nat- 
urally explained. Such inhomogeneities are taken into 
account by a charge distribution, but as discussed above, 
the main features of our calculations are independent of 
the details of the probability charge distribution. In con- 
clusion it is crucial to take into consideration the fact 
that the HTSC are inhomogeneous materials in order to 
describe the main qualitative features of the high-T c su- 
perconductors. 
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